function [c,ceq] = outputcon(sbar,eq,p)

stilde = eq(1);
xlow = eq(2);
xup = eq(3);

r = laborassignment(sbar,stilde,xlow,xup,p);

c = p.contol + eps - r.c; % ensures that c violates c<=0 if r.c=0 and p.contol is the tolerance for constraint violation
ceq = [];

end